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Abstract: In this article we investigate charged particles in gauge theo- 
ries. After reviewing the physical and theoretical problems, a method to 
construct charged particles is presented. Explicit solutions are found in 
the Abelian theory and a physical interpretation is given. These solutions 
and our interpretation of these variables as the true degrees of freedom 
for charged particles, are then tested in the perturbative domain and are 
demonstrated to yield infra-red finite, on-shell Green's functions at all or- 
ders of perturbation theory. The extension to collinear divergences is stud- 
ied and it is shown that this method applies to the case of massless charged 
particles. The application of these constructions to the charged sectors of 
the standard model is reviewed and we conclude with a discussion of the 
successes achieved so far in this programme and a list of open questions. 
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1 Introduction 

One of the most significant advances in particle physics was the definition of what a 
particle should be. Wigner |IJ identified a particle with an irreducible representation 
of the Poincare group or its covering. Such representations have a well defined mass 
and spin. Building upon free creation and annihilation operators, and invoking |2| the 



1 



cluster decomposition theorem, one can construct the paradigm quantum field theoretic 
description of particles scattering into other particles. 

Physical particles generally carry additional quantum numbers such as isospin, elec- 
tric and colour charge. These particles are described by the gauge theories of the 
standard model. Due to the masslessness of the gauge bosons, these field theories 
are plagued by infra-red divergences which dramatically alter the singularities of the 
Green's functions of the matter fields. The states no longer form irreducible represen- 
tations of the Poincare group and it is therefore widely believed |4[] that there is 
no particle description of charges such as the electron. This review will describe the 
physics underlying this problem, the manifold consequences of this breakdown and how 
a particle description of charged fields may be recovered. 

The standard picture of a particle's journey to a detector is as follows. When it is 
shot out of a scattering event it is, after enough time, a long way away from any other 
particle and the interaction between it and the other particles may be neglected. The 
free Hamiltonian then describes its dynamics and a particle description holds. This, 
the cornerstone of the interaction picture, breaks down for charged particles interacting 
via the gauge theories of the standard model. 

The importance of the large distance interactions, which mean that the residual 
interactions between particles at large separations cannot be neglected, is most obvi- 
ous in Quantum Chromodynamics (QCD). Experimentalists observe hadrons and not 
quarks which are professed to be eternally confined. Although hadrons are colourless 
(chargeless as far as the strong nuclear force is concerned) their masses are rather 
well described in terms of building blocks which are the original Gell-Mann quarks. 
However, the current quarks of the QCD Lagrangian are not the constituent quarks of 
hadronic spectroscopy. This dichotomy is an unsolved puzzle and shows itself in many 
ways. For the light quark flavours the masses of the current and constituent quarks are 
very different (roughly two orders of magnitude for the u and d flavours). The division 
of the spin of the proton amongst its constituents has led to the 'proton spin crisis', 
which name shows the depth of our difficulties in understanding the experimental data. 
Finally, since constituent quarks are presumably constructs made from surrounding a 
Lagrangian matter field with a cloud of coloured glue, it is initially at least highly 
unclear how the constituents obtain a well-defined colour. 

A major aim of the programme of research |J described here is to understand how 
constituent quarks arise in QCD, their mutual interactions and finally how it is that 
these effective degrees of freedom are confined. 

Electrons are not confined and so it may seem initially less obvious that the inter- 
action picture paradigm breaks down here too. Yet the S-matrix elements of Quan- 
tum Electrodynamics (QED) are afflicted by infra-red divergences, as are the on-shell 
Green's functions of the theory with external legs corresponding to charged particles 
(which we will sometimes generically refer to as electrons). The infra-red problems 
change the form of the Green's functions such that we cannot associate a pole to the 
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external legs. There are two different sorts of divergence here: soft divergences, which 
show up in Green's functions and S-matrix elements, and phase divergences which 
occur in the phase. These latter divergences are often ignored in QED, but are of 
importance in QCD (see Sect. 3.4 of Ref. || and also 0). 

The masslessness of the photon is the underlying cause of the infra-red problem^. 
This vanishing of the mass means that photons can travel over a large distance and 
indeed that an infinite number of soft photons can be created for any finite amount 
of energy. Recognition of this led to the Bloch-Nordsieck answer to the infra-red 
problem in QED: since any finite experimental resolution does not restrict the number 
of photons which may accompany any charged particle, an experimental cross-section 
must come from summing over all these possibilities. In this sense QED is taken to 
be a theory defined only at the level of (measurable) cross-sections and not in terms 
of (unobservable) S-matrix elements. Although this is in no way incompatible with 
experiment, it is a radical conclusion. 

The survival of large distance interactions is responsible for the claim |3J that 
charged particles cannot be incorporated into relativistic quantum field theory. This 
conclusion followed from noting that as the usual free Hamiltonian does not determine 
the asymptotic dynamics, it must be modified. This leads to a description in terms of 
coherent states and there is no pole structure. Is particle physics then a misnomer? 

In fact we will see below that these problems can be resolved by realising that, as 
the asymptotic interaction Hamiltonian is not zero, electric charges are surrounded by 
an electromagnetic cloud, just as quarks are by glue. It will then be seen that such 
systems of dressed charges must be gauge invariant. We will demonstrate that the 
Green's functions of such dressed charges are infra-red finite and have a good pole 
structure. 



1.0.1 General Properties of Charges 

The gauge dependence of the Lagrangian degrees of freedom means that it is hard 
to associate any physical meaning to them. This hinders any attempt to understand 
how phenomenological models and concepts can arise from the underlying theory. As 
an example of this we recall that the highly non-trivial gauge dependence of the La- 
grangian vector potential in non-abelian gauge theories has put difficulties in the way of 
extracting such, phenomenologically useful, ideas as effective gluon masses from lattice 
calculations || . Constrained dynamics is the mathematical tool appropriate to finding 
gauge invariant degrees of freedom. The true degrees of freedom then correspond to 
locally gauge invariant constructs (which obey Gauss' law). There have been many at- 
tempts to obtain such variables in gauge theories (see, e.g., and references therein). 

4 For this reason one sometimes speaks of mass singularities. 
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We now want to sketch some general properties that any description of charges must 
fulfill. 



Local gauge transformations in QED have the form 

A^x) -> A„(x) + d^x) , and $(x) -> e ie6{x] } ip(x) , (1) 

so if the coupling e could be switched off, then the Lagrangian fermion would be locally 
gauge invariant. Similarly Gauss' law 

d l F l0 = -eJ , (2) 

where Jq is the charge density, would, in the e = limit, imply that only the transverse 
components of the field strength were physical. However, as we have noted the coupling 
does not vanish and so the matter fields cannot be identified with physical particles. 
Similarly Gauss' law shows that there is an intimate link between the matter fields and 
the electromagnetic cloud which surrounds them. An immediate consequence of this 
is that any description of a charged particle cannot be local since the total charge can 
be written surface integral at infinity. 

We also see that objects which are invariant under global gauge transformations 
are chargeless, as the charge density is the generator of such global transformations. 
The Gauss law constraint generates local gauge transformations, so we demand invari- 
ance under local but not global gauge transformations of any description of a charged 
particle. 

The form of the cloud around a charge determines the electric and magnetic fields 
surrounding the charge. This implies a fundamental non-covariance: the velocity of 
any charged particle will determine the nature of the cloud. 



The implications of these inevitable properties of non-locality [T(J and non-covariance 



TT|| for any description of charged particles for the general properties of gauge theo- 



ries have been investigated both for scattering theory |T2|, [13j and in axiomatic ap 



proaches 14, 15 



1.0.2 Dirac's Dressings 



To the best of our knowledge the interplay between a charged matter field and the 
electromagnetic cloud which inevitably surrounds it was first used by Dirac to 
construct a specific description of the electron. He suggested that one should use 



i) D {x) = exp -ze-^ ip(x) . (3) 
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This he motivated in the following way: it is locally gauge invariant^. Using the funda- 
mental equal-time commutator, [Ei(x),Aj(y)] = i5ij5(x — y), and the representation 

d i A i, s 1 f,3 diAi(x ,y) 

-=z-(x) = ay — ; ; — 4 

V 2 V ' 4vr J y \x-y\ v ' 

then the electric field of the state ipu(x)\0) is found to be 

E\x ,y)Mx)\0) = ~ ^%^)|0>, (5) 

47r \x — y \ 6 

which is what one would expect for a static charge. He further pointed out that this 
is actually a member of an entire class of composite fields 

= exp (-ie J d 4 zr(x - z)A fl {z)^j ij>(x) , (6) 

which, he argued, are gauge invariant for all if it is demanded that d^f^iw) = 5^(w) 
holds. These constructs are all evidently non-local. We refer to the cloud around the 
matter field as a dressing. It is clear that the dressing suggested by Dirac (|3]) is 
also non-covariant. This review will be concerned with generalising and refining such 
descriptions of charges. 

Such pictures have been rediscovered by various authors since Dirac and there have 
been many attempts to use certain examples of this wide class of dressings over the 
years |T7|, [0§ 0, |21|, |22|, [2|, JHj, |2|, g§ 0, [§, |2j| |3g . 



A natural extension of this is to generalise Eq. |^ to systems involving more than one 
matter field. Two opposite charges at different points can be made gauge invariant by 
including a dressing which keeps the entire system gauge invariant. Such a description 
could correspond to a positronium state, a hydrogen atom or a meson. If the cloud 
factorises so that each of the matter fields together with its part of the dressing is 
gauge invariant, then we can clearly speak of constituents. A dressing which does 
not factorise at all would mean that we had an effective 'meson' field but could not 
really speak of constituent particles. The success of the constituent quark picture of 
hadrons, taken together with the rough equality of the constituent masses in mesons 
and baryons indicates that there must be at least a rough factorisation of the dressings 
in some dynamical domains of QCD. 

Once one has an ansatz for the dressing, it may be used to study the interaction 
between charges. The easiest manner to produce a gauge invariant description is to 
link the two matter fields by a path ordered exponential along some line. The potential 
is then obtained^ by taking the commutator of the Hamiltonian with this description: 
that part of the energy which depends on the separation is the potential. The string 



5 But not globally gauge invariant. 
6 For details, see Sect. 7 of § and 
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ansatz leads already in QED to a linear, confining potential. Furthermore the overall 
coefficient 1S ; cLS cL result of the infinite thinness of the string, divergent. This string 
description has no physical relevance, even in QCD the potential between two heavy 
quarks should have a Coulombic form at short distances [M] and the finite string 
tension implies a cigar shaped dressing. 



Since the string model corresponds to an (infinitely) excited state, it is unstable PT 



If we consider two extremely heavy, fixed charges and neglect pair creation, we may use 
the free Hamiltonian and thus solve the time development of the system starting from 
the string ansatz initial state. The electric and magnetic fields immediately broaden 
out and, of course, lead to the usual Coulombic field. An animation and detailed 
discussion of this can be found at the web site: http: //www. if ae . es/~roy/qed.html 

After this fly- by tour of the subject, it is time to get to grips with the physics of 
charges. The structure of this article is as follows. In Sect. 2 the form of the interaction 
Hamiltonian at large times is investigated. It is shown not to vanish for the matter fields 
but we demonstrate that, for suitably dressed fields, the asymptotic interaction does 
indeed vanish. This result is then used to construct explicit charged fields. In Sect. 3 
an alternative account of the charged fields in the heavy scalar theory is presented. 
These results and our interpretation of them are then tested in perturbation theory in 
Sect. 4. It is shown that the soft and phase divergences cancel and a pole structure is 
obtained. The massless electron limit is considered in Sect. 5 as a testing ground for 
the study of collinear divergences. The form of the asymptotic interaction Hamiltonian 
is shown to be such that the solutions of the dressing equation will remove these mass 
singularities as well. Subtleties of the dressing equation in this limit are investigated 
in this section. Finally Sect. 6 reviews what we have learned and presents a list of 
outstanding problems. 



2 Charged Particles 



As we have seen, the masslessness of the photon implies that the interactions between 
charges cannot be 'switched off' in the remote past or future. It is this long range 
nature of electromagnetic interactions that lies at the heart of the infra-red problem 
in QED. In this section we will review how the asymptotic dynamics found in such 
a gauge theory deviates from that of a free theory. As a consequence of this we will 
see that the correct identification of asymptotic, charged particle states can be made 
through a process of dressing the matter of the theory. 



2.1 Asymptotic Dynamics 

The S-matrix codifies the intuitive picture of a scattering experiment whereby in- 
coming particles get sufficiently close to interact, resulting (after the dust settles) in 
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some set of out-going particles. The in-coming and out-going particle regimes are iden- 
tified with elements of the Fock space constructed out of the creation and annihilation 
operators found in the free theory. As discussed in the introduction, the reason for 
this is that it is these states that can be identified with (tensor products) of irreducible 
representations of the Poincare group. As such, they are particles! 

The key assumption, then, in the S-matrix description of particle scattering is 
that in- and out-regimes where the dynamics is that of the free system exist, i.e., the 
interacting Hamiltonian must tend to zero in the remote past and future. It is this 
assumption that fails in a gauge theory |32|. || . 



To see how this arises in QED, we start with the gauge fixed Lagrangian density: 
£qed = -f^ + #7 M («9 M - ieA^ - + \^B 2 + d^B . (7) 

In this we are working in the Lorentz class of gauges (Feynman gauge corresponding 
to the choice £ = 1). The gauge invariance of the physical states is encoded in the 
condition that _B + |phys) = 0. For the present, we do not allow the mass of the electron 
to be zero (m ^ 0). 

The interaction Hamiltonian is given by 

H hit (t) = -ej d 3 x A^t, x)J"(t, x) , (8) 

where the conserved matter current is J M (t, x) = ip(t, x) r j IJ "ijj(t, x). In order to construct 
the S-matrix, the fields that enter this part of the Hamiltonian are taken to be in the 
interaction picture. We recall that this means that the time evolution of the states is 
described by (H) while the fields evolve under the free Hamiltonian. Thus in @ we 
should insert the free field expansions for the matter and gauge fields, this will then 
allow us to study in detail the large t behaviour of the interaction. 

We take as our free field expansions 

= "y= {bfr, s)u s (p)e-^ + dt(p, S )t;'(p)e**} (9) 



and 



M*) = J — {a M (fe)e-*« + aV(fc)e^} , (10) 



where E p = y\p\ 2 + m? and Uk = \k\. Inserting these into (|8|) results in eight terms 
which we group according to the positive and negative frequency components of the 
fields. Each of these pieces will have a time dependence of the form e tat where a involves 
sums and differences of energy terms. As t tends to plus or minus infinity, only terms 
with a tending to zero can survive and thus contribute to the asymptotic Hamiltonian. 
After performing the spatial integration, and using the resulting momentum delta 
function, only terms of the form e ±tt ^ E p+ k ~ Ep±L ° k ' > have a large t-limit; there are four of 
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them. The requirement that E p+ k — E p ± ujk ~ can only be met in QED because the 
photon is massless, in which case it implies that Uk ~ 0, i.e., only the infra-red regime 
contributes to the asymptotic dynamics. From this observation it is straightforward to 
see that the full interacting Hamiltonian (H) has the same asymptotic limit] 7 ] as 

7C(t) = - e / d " x M*> X ) J ^ x ) ( n ) 

with 

d 3 p p^ _ / ,_^ 3 / p 



■^'•-Wc^i^n-tv- (12) 

The operator content of this current is only contained in the charge density 

p(p) = J2( b ^(p> s ) b (p> s ) ~ rft (^' s ) d (p> s )) ( 13 ) 

s 

which implied that the asymptotic current satisfies the trivial space-time commutator 
relation 

-C(v)] = o. (14) 

As such, this current can be interpreted as effectively the integral over all momenta 
of the current associated with a charged particle moving with velocity p^/Ep. Such a 
current does not vanish as t — > oo. We thus see that the asymptotic dynamics of QED 
is not that of a free theory. 

The non-triviality of the asymptotic dynamics dramatically alters the form of the 
in- and out- matter states. In particular, their propagator no longer has a pole-like 
structure, but instead behaves near its mass-shell like (p 2 — m 2 ) 13 " 1 where the exponent 
(3 is gauge dependent; in the Lorentz class of gauges it is given by 

£ = 3). (15) 

It is this observation that lies at the heart of the statement [§J] that there is no relativistic 
concept of a charged particle. 



2.2 Charges as Dressed Matter 

The persistence of the asymptotic dynamics in QED means that we cannot set the 
electromagnetic coupling to zero for the in-coming and out-going particle. As an im- 
mediate consequence of this we see that the matter field, ip(x), cannot be viewed as the 
field which creates or annihilates charges since it is not gauge invariant in the remote 
past or future. An equivalent statement of this fact is that the matter field, ip(x), does 
not satisfy Gauss' law at any time. 

7 A fuller account of this derivation of the asymptotic interaction will be presented in Sect. 5 as 
part of a more general analysis that includes the case of massless matter. 
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To construct a charged field we need to be able to find a functional of the fields, 
h~ 1 (x), such that, under a gauge transformation described by the group element U(x) = 

e ieO(x) ^ we J^yg 

h~ 1 {x) -»• hr\x)u(x) . (16) 

Then, from (JJ) the charged field will be given by the gauge invariant product 

= h~ 1 {x)if){x) . (17) 

It is this product which makes precise how we associate charges with dressed matter: 
a (chargeless) functional of the fields which transforms as in ( |T6|) is what we mean by 
a dressing; the product ( |TTD we then identify as charged matter. We continue to refer 
to the matter terms that enter directly into the Lagrangian simply as matter. 

There are many distinct ways to construct dressings, reflecting the fact that our 
identification (|T7| ) of charged matter should be viewed as a minimal requirement. In 
specific applications the form of the dressing must be tailored to the physics at hand: 
charges that enter into bound states will have a very different structure to those that 
describe particle. In our programme to construct charges explicitly we have so far dealt 
exclusively with the construction of charged particles^]: it is an open and immensely 
interesting problem to extend our approach to bound states. 

In order for the charged matter to describe a particle, we demand that the dressing 
is such that the charged field creates a state which is an eigenstate of the energy- 
momentum tensor, i.e., that it has a sharp momentum. For the matter field, if), this 
was not possible due to the infra-red problem ||. For our charged field, though, we 
shall see that the form of the dressing can be chosen so that there are no infra-red 
divergences. 

The interacting Hamiltonian (H) is derived from the coupling term in the matter 
part of the Lagrangian density for QED: 

#(x) 7 M (^ - ieA„(x))i/>{x) . (18) 

We can rewrite this in terms of the physical charged fields *f>(x) as 

itf (aOy(d„ - ieAl(x))V(x) (19) 

where 

Al = h- l A»h + U> il {h- 1 )h, (20) 

which we recognise as a (field dependent) gauge transformation?] of the vector potential. 
Thus, written in terms of the charged fields, the interacting Hamiltonian is 

HL(t) = -e / d 3 x A% x)J^t, x) , (21) 

8 In non-abclian gauge theories the gauge field also needs to be dressed in order for it to carry 
colour HHH. 

9 The non-abelian appearance of this gauge transformation in QED arises from the fact that the 
dressing is field dependent and thus might not commute with the potential. 
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which asymptotically becomes 



HL(t) - -eJd 3 xA h ^t,x)J^t,x) (22) 

r d 3 p A h ^(t,x)p^ 3 / p \ 
dX Ij^f ~^^ P{P)6 [ X ~ E?) ■ (23) 

This would vanish if we could construct the dressing such that A^(t, x)p fl = 0. 

To investigate the extent to which this asymptotic interaction can vanish, we note 
that the momentum, is an arbitrary on-shell, four vector. Now through a gauge 
transformation, such an algebraic condition can be imposed on the vector potential at 
one point in the mass shell, but not for the whole mass shell. So we cannot expect to 



be able to construct a dressing such that the asymptotic interaction Hamiltonian (|22|) 
vanishes for the whole mass-shell. However, we can insist that our charged field creates 
an in-coming or out-going particle with a definite momentum. 



From this discussion we see that if we want the charged field fllTD to asymptotically 
create or annihilate a charged particle moving with four-velocity then the dressing 
must satisfy the additional kinematical condition that 

t*Mj(f, a) = . (24) 

This will then ensure that, at the point in the mass-shell where p^ = mu^, the asymp- 
totic interaction Hamiltonian (|22"D vanishes and thus the state created by the field will 
have the appropriate sharp momentum. 

What this means in practical terms is that each charged particle must be con- 
structed out of the matter fields with a different dressing — reflecting the velocity of 
the particle concerned. At first sight this looks very peculiar: we are used to particles 
being put on-shell as the mass of a particle is a well defined quantum number, but a 
particle's velocity is not usually thought of as a quantum number of the system. For 
charged particles, though, velocity is a well defined quantum number. In order not to 
interrupt our account of how to construct charges we shall postpone a discussion of 
this interesting point until Sect. 3 where the connection with the heavy charge sector 
will also be discussed. 



2.3 Construction of Charges 

In order to construct an asymptotic charged particle moving with four- velocity u^, we 
need to be able to find a dressing field h~ x such that flT6|) and (p4] ) hold. The dressing 
equation (|24]) is, for this purpose, best written as 

(77 + vYd^hT 1 ) = -iehT 1 ^ + vYA^ , (25) 
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where we have written the four- velocity as 7(77 + t>) M , where 77^ is the time like vector 
(1, 0), v M is the space-like vector (0, v) with v the three- velocity of the charged particle 
we wish to construct and 7 is just the standard relativistic factor (1 — v 2 )~z . 

To see how to construct the dressing, we first restrict ourselves to lowest order in 
the coupling. That is, we take = 1 — ieR(x) and, to this order, the dressing 

equation is 

(77 + vydpRix) = (r] + vfA^x) . (26) 

This equation can now be easily solved in terms of an integral along the world line of 
a particle moving with the given velocity: 

R{x)= {v + vfA^dt + xixa)- (27) 

J a 

In this expression we have, for convenience, parametrized the world line not by the 
proper time but by xf = x^ + (t — x°)(r] + v ) fl and we have introduced an arbitrary 
reference time a. The term x( x a) is> f° r the moment, an unspecified field configuration 
in the kernel of the differential operator (77 + v)^d^. It is useful here to group the 
a-dependent terms together under the integral, and write R(x) as 

R(x) = £ (( v + vyA^x t )-^(x t )yt + X (x) 

= [ (v + vr(A^x t )--^(x t )yt + X (x). (28) 

In addition to the dressing equation, we also need to ensure the correct gauge trans- 
formation properties for the dressing. To this order in the coupling, this means that 
under the gauge transformation A^ — > A^ + d^d, we must have 

R(x) -> R(x) + 0(x) . (29) 

This, in turn, implies that 

X (x) -> X (x) + 9{x) . (30) 

This is precisely the type of transformation that Dirac investigated in (|6]) and we would 
be tempted to write 

X (x)= Jd i zr(x-z)A,(z), (31) 

where f^{x — z) satisfies d il f >1 {x — z) = 5 4 (x — z). But this only implies ( |3T)D if no 
surface terms arise when we integrate by parts. The restriction on the local gauge 
transformations to those that vanish at spatial infinity is quite natural as finite energy 
restrictions impose a 1/r fall-off on the potential. However, no such restriction applies 
to the fields at temporal infinity. Thus the form of f^{x — z) must be such that no 
such surface terms arise. 
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As it stands, we can only infer from this that f°(x — z) should be zero outside of 
some bounded region in the z°-direction. To get more from this, we make the ansatz 
that in fact 

X (x) = J d 4 z G(x - z)G^A ll (z) , (32) 

where is a first order differential operator and Q ■ dG(x — z) = 5(x — z). In order 
to avoid the surface terms that would obstruct the gauge transformation properties of 
the dressing, we must have that the operator Q ■ d cannot involve any time derivatives. 
Given this restriction, we see that G(x — z) = 5(x° — z°)G(x — z). We shall, for 
convenience, write \ as 

g . j\ 

The temporal parameter a which enters into the form of the dressing fl2~7|) has been 
introduced by hand and thus should not directly affect any physical results. Taking 
the derivative of ([28|) with respect to a yields the gauge invariant term 

- (rj + vr(A,(x a ) - dtfM) = -(V + vY(A„(x a ) - ■ (34) 

This we can write in a form where the gauge invariance is manifest as 



-tn + vr[^f^(x a ) - ) = -(v + vY^M . (35) 

The condition selecting physical states is that they are annihilated by the 5-field, thus 
any physical observable must commute with B(x). The relevant equation of motion 
that follows from ([?[) is, to this order, d u F Ufl = d^B. This tells us that, in order for 
(35) to act trivially on physical states, it must be equal to 

(v + vT^(x a ). (36) 

This now allows us to find the form for and hence the dressing to this order. 

The first order operator is constructed out of the vectors <9 M , n M and that 
characterise the theory. The anti-symmetry of the field strength F UfJ- in (|35| ) implies 
that (RBI) will follow if 



g» = -d v + (77 + v) u (a(n-d) + P(v-d)) . (37) 

In which case the second order operator Q ■ d is 

g . d = V 2 + (3{vd) 2 + {a- l){n-d) 2 + {a + /3)(v , (38) 

where V = didi. We have seen that it is essential for the gauge invariance of the charges 
that this operator has no time derivatives in it. This follows only if a = 1 and (3 = — 1. 
Hence we see that 

g» = { rj + v y{n-v)-d-d» . (39) 
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Thus, to lowest order in the coupling, we have seen that h (x) = 1 — ieR(x) where 
R(x) is the sum of two terms: 

r x° qh p G^A 

R{x) = -j a ds( V + vT-^f{x t ) + Y^ {X) - (40) 

As will become apparent when we perform perturbative test of our construction, this 
decomposition of the dressing into two terms reflects the two manifestations of the 
infra-red in QED that were discussed in the introduction: the soft divergence and the 
phase divergence. Indeed, the first (gauge invariant) part of R(x) will be shown to 
be responsible for controlling the phase structure of charged matter; while the second 
term will be responsible for the soft dynamics of the charged particle. 

It is, perhaps, helpful to specialise to the static situation in order to get a better 
feel for the form of the dressing we have been constructing. In that case v = and we 

get v( ^ r\ dip ^t ^ d ^t ^ 

R{x) = - ds^-(s,x) + —(x), (41) 

where 1/V 2 has been defined in (|4]). The expected spatial non-locality of the charge is 
now manifest. More generally, the inverse to Q ■ d is given by 

'/M^-ih/H- (42) 



Q ■ d An J \\x — y 

where 



d'k — -r , (43) 



2vr 2 J V v ■ k 



and = (t) + v)^{rj — v) ■ k — 



In order to go beyond this first order result we must take into account the fact 
that the fields which enter into the dressing are operators, and thus we need to know 
their commutation properties. Equal time commutators between the fields can be 
directly read off from the Lagrangian (0). However, the non- locality displayed in (|40[) 
means that we need general space-time commutators between the fields. For the fully 
interacting theory there is no general approach to finding these commutators without 
first solving the whole theory. Charges as particles, though, only make physical sense 
in the asymptotic regime governed by the interacting Hamiltonian (pi] ). In this regime, 
the relevant dynamics of the gauge field is described by the simpler set of equations: 

d v F Vil = d^B-eJ™ (44) 
^ = -£B, (45) 

where the matter coupling is now through the asymptotic current ( |T2"| ) . For this theory 
the space-time commutators can be explicitly constructed, as we will now show. 
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In Feynman gauge (£ = 1) the potential satisfies 



DA, = -eJ™ . (46) 

The general solution to this is 

A,(x) = A^(x) - e J d*yD R (x - y)J™{y) , (47) 

where Dr(x — y) is the retarded Green's function for the □ operator and A f * ee is a 
solution of the homogeneous equation DAjf 26 = 0. From the form of this solution for 
the potential, and given the triviality of the asymptotic current commutations (|14|), 
the space-time commutator for the potential in this interacting theory are seen to be 
the same as that for the free theory. 

The free field equation n^4r ee = implies that 



fr ee ( x ) = f d 3 z{d z D{x - z)Af e (z) - D(x - z)d Q Al cc (z)) (48) 



.4 

where D(x — z) is the commutator function for free fields: 

/d^k 1 
j^r ^f Mx ~ v) sin ( Wfc(x ° - ^ 0) ) • (49) 

The identification in fl4"8] ) is made by first observing that the right hand side is indepen- 
dent of z°: setting z° = x° then implies the result. Exploiting this ^-independence, the 
commutator [A^ rce (x), ^4^J ee (2/)] is simply calculated by using (|^) with z° = y°. Then 
the equal time commutation relations [A^y), A v (z)] et = —ig, u 5(y — z) can be used. 
This, in conjunction with our observation that the free and asymptotic interacting 
fields have the same commutators, results in the space time commutators in Feynman 
gauge being: 

[A,(x), A u (y)] = [A^(x), A^y)} = -ig»„D{x - y) . (50) 

Using this, the first order result (f40"D can be exponentiated (while still preserving the 
a-independence of the construction on physical states) to give 

h-\x) = e- ieK{x) e- iex(x) , (51) 

where 

x(*) = fr7f(*). (52) 

and 

K(x) = -\ a ds(r, + vy^f{x s )+e j_Js( V + v),—(x s ) (53) 

" e7 " u ' x ' J f (2^ e - ( k . v y (54) 
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These results are, by construction, independent of a. Setting a = — oo gives (modulo 
the field independent, tadpole term) 



K(x) = - ds ( V + v r^f(x s ) , (55) 

which is a form adapted to performing perturbative calculations. The full derivation 
of these results will be presented elsewhere. 

In summary, we have seen that charged fields corresponding to charged particles 
moving with momentum, p M = mw M may be described by 

* p (x) = e- teK " {x) e- tex " {x) ^(x) , (56) 

where the dressing satisfies both the gauge transformation property ( |T6"D and the dress- 
ing equation ( p5|) with the four velocity appropriate for the particular momentum p. 
Our construction has shown that the on-shell Green's functions of these dressed fields, 
taken at the correct points on their mass shells, should not suffer from any infra-red 
divergences. Such formal arguments must be tested: in Sect. 4 we will show that this 
prediction || is true to all orders in perturbation theory |36] . 



We now conclude this construction of charges in QED with an observation and a 
comment. Note that in the static limit the charge field is given by 

exp Ue y ds—^-{s,x) J exp I — ze— ^- I ip{x) = exp \ ie J ds-^-{s,x) I ip D {x) 

(57) 

where ipo was Dirac's proposal @ for a static charged field based on the form of the 
electric field produced by such a charge. The additional term in (|57D, that follows 
from the more fundamental dressing equation approach presented here, does not alter 
the form of the electromagnetic field and was thus missed by Dirac. Its role in the 
cancellation of the phase divergence will be presented in Sect. 4. 

The construction of the charge above has been done in Feynman gauge. It is es- 
sential to verify that the results hold for the full Lorentz class of gauges. While it is 
true that the very construction of the charges (|i~7|) ensures gauge invariance, it is not 
immediately clear that this is equivalent to the independence of the construction from 
the gauge fixing parameter £. In particular, changing £ from one will modify the form 
of the space-time commutators, fl5PD, and thus could have a potentially non-trivial 
impact on the construction of the dressing. A full proof of the gauge invariance of 
these dressed fields will be presented elsewhere, however, in Sect 4, we shall demon- 
strate perturbatively that the n-point Green's functions constructed out of the charged 
fields are manifestly ^-independent. Before making these tests, we will now present an 
alternative approach to the physical restrictions on dressings. 
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3 Heavy Charges 



We have seen in the last section that the construction of the dressing for a charged 
particle depends on the velocity of the charge. In this section we wish to discuss, in 
more general terms, the role the velocity of a charged particle plays in the structure of 
the states of QED . We will then look at heavy charges, where velocity is also singled 
out, and see how the asymptotic regime, governed by the asymptotic Hamiltonian (PD , 
can be characterised in terms of the mass of the charge. In addition, we will also see a 
further derivation for the dressing equation (f25|). 



3.1 Moving Charges 



It is not possible J57J] to talk about the space of states of QED without first identify- 
ing the observables of the theory and their algebra. The dual requirements of gauge 
invariance and locality on any physical observable have immediate and quite striking 
implications for the structure of QED. 

A familiar example of this is the fact that charge is superselected: the state space 
is a direct sum of different charge sectors and the action of any physical observable 
cannot change the charge. In order to derive this physically important result we recall 
that Gauss' law (0) implies that on physical (gauge invariant) states there is a precise 
relation between the electric field and the charge density. In particular, the total charge 
Q is given on such states by the spatial integral 



Q= fd 3 xdiEi(t,x). (5£ 



Using Gauss' theorem we can relate this expression for the charge to a surface integral 
and, in particular, to one at spatial infinity. That is, if we define the electric flux in 
the direction x by 

£i(x) := lim R 2 Ei(x + Rx) , (59) 

R—>oo 

then the total charge is given by the 'flux at infinity' through 



Q= ds-S. (60) 

Js*, 

The superselection properties of the charge then follow immediately from the obvious 
fact that any local observable will always commute with the flux at spatial infinity. 

This simple argument can be extended to show why velocity plays such a central 
role in the identification of charged particles. Recall that the the electric field of a 
charge whose present position is y and moving with velocity v is 

7 2 [(v ■ (x - y)\ +\x-y\ 2 
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The electric flux at spatial infinity is then 



Just as above for the total charge, this flux will commute with all local observables 
and thus can be used to characterise distinct sectors of the theory. As this flux only 
depends on the velocity of the charge we see that, in the asymptotic regime where 
charged particles exist, their velocity is a well defined quantum number. 

Within the class of local observables, this means that velocity is superselected. 
However, the very non-locality of the construction of the charges implies that we have 
to extend the algebra of observables to include non-local constructions such as dressings. 
As we will see in the next section, the perturbative construction of Green's functions 
for charged matter will involve non-local interactions between different charges and 
thus allow charges with different velocities to interact. This does not mean, though, 
that the total charge Q is no-longer superselected: the dressings are all chargeless and 
hence the interactions they induce commute with the charge. 



3.2 Heavy Matter 

In our discussion of the asymptotic interactions found in QED, we have simply taken the 
large t-limit of the interaction Hamiltonian in the interaction picture. A disadvantage 
of studying this limit, though, is that it does not supply a description of the scale for 
the onset of the asymptotic dynamics and thus of the domain within which velocity 
is a valid quantum number. More properly, then, we should investigate the limit as 
some dimensionless parameter gets large. In QED with massive matter, the natural 
parameter is the product of t and the mass scale of the theory set by m, the lightest 
mass of the system. In this sense, the dynamics at large time is equivalent to that of a 
theory with heavy charges. Given the central role of the asymptotic dynamics in our 
programme, it is important to understand precisely how the asymptotic interactions 
emerge in the heavy sector, and to understand the significance of the dressing equation 
for heavy matter. 

It is well known || that the infra-red structures found in QED are independent 
of the spin of the charged particles, thus it is instructive to see how the interaction 
(|TT| ) also emerges from the heavy sector of scalar QED. The matter part of the QED 



Lagrangian is now 

{Pp<f>y(p»<i>)-m 2 <ft<i>. (63) 



The heavy limit can only be taken at specific points on the mass-shell of the 



particle. To this end, one introduces the rescaled fields 

4>{x) := V2^e imu - X <f){x) , (64) 



17 



where we have chosen a four-velocity (u ■ u = 1) that will ultimately describe the 
velocity of the heavy charge. In terms of these new fields, the matter part of the 
Lagrangian becomes 

iftvTD^ + ^-{D^)\D^) . (65) 

In the large m-limit only the first term survives and the equations of motion for the 
heavy matter become 

u»D^ = . (66) 

The interaction Hamiltonian is easily identified in this limit and is constructed out of 
the current 

^(1)=^!), (67) 

which has, as expected, precisely the same form as the asymptotic current (|12|) for the 
specific point on the mass-shell described by . 

Heavy charged matter also needs to be constructed out of the heavy matter through 
the process of dressing. Thus the heavy scalar charge is given by the gauge invariant 
field 

= h' l {x)4>{x) , (68) 
where the dressing is that appropriate to a charge moving with the given four velocity 

vPdpihr 1 ) = -ieh-^Ap . (69) 

This dressing equation, in conjunction with the equation of motion for the heavy matter 
(|66| ) implies that 

uUdfi = . (70) 

That is, the state created by the field (f = e~ imu ' x $> is an eigenstate of the momentum 
operator: 

i>|0)=m W |0), (71) 
i.e., as claimed, it is a particle with momentum mw M . 

Having seen that the heavy sector characterises the regime described by the asymp- 
totic Hamiltonian (PD, we can now identify the onset of the asymptotic dynamics, and 
thus the emergence of a charged particle interpretation, with the kinematical regime 
where the lightest mass m of the matter dominates the momentum flow in any process. 



This is precisely the domain of the eikonal approximation in soft photon amplitudes [39 



whereby the photon has a spin-independent coupling to the matter, which acts as an 
external point-like source. It is also the domain where scattering theory applies and 
we will now test our construction. 
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4 Perturbation Theory 



Our aim now is to show how we may apply the physical, dressed fields we derived in 
the last section to remove infra-red divergences already at the level of Green's func- 
tions. We will start with a one-loop example, where we will classify the various sorts of 
mass singularities present in an unbroken Abelian gauge theory. Then we will demon- 
strate |36| their explicit cancellation when physical, dressed fields are used. Finally we 
will show that this cancellation is in fact general and occurs at all orders. In this review 
we will solely concern ourselves with the IR-structure and not with the renormalisation 
of UV-singularities — full calculations of the renormalisation constants associated with 



the physical propagators of spinor and scalar QED were, however, given in Ref.'s 40 



and [fy]] respectively. Ref.'s P, |42|, ^9f offer introductions to infra-red divergences in 
perturbation theory. 



4.1 One-Loop Calculations 

The full spectrum of infra-red divergences that characterise gauge theories where the 
charged particles are massive all occur in pair creation. We thus start by briefly 
reviewing this process for undressed matter fields. We will use the theory of scalar 
QED since these divergences are insensitive to the spin of the particles and we may so 
avoid the technical complications of Diracology. The Feynman rules are 




Fig. 4.1: The Feynman rules. 

We have chosen to work in Feynman gauge, but the physical results will be shown 
to be gauge invariant below. The seagull (<ft*<f)AA) vertex does not alter the spin 
independent IR-behaviour and is therefore irrelevant to our purpose. 



4.1.1 One-Loop Pair Creation 

Consider now pair creation from a classical source^] as shown in Fig. 4.2. 

10 The Feynman rule for the vertex with the source is here taken to be unity for simplicity, it could 
also be given the renormalisation group invariant form, to 2 </>*0. 
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Fig. 4.2: Covariant pair production diagram. 

The vertex function corresponding to Fig. 4.2 is given in Feynman gauge by 
-ie 2 f d 4 k (2q + k)^(2p - k) v g„ v 



[p 2 — m 2 ] [q 2 — m 2 ] J (27r) 4 [(q + k) 2 — m 2 + ie] [(p — k) 2 — m 2 + ie] k 2 + ie 



• (72) 



Naive power counting tells us that this does not have an infra-red divergence when the 
outgoing particles are off-shell. However, if we extract a simple pole for each of our 
outgoing particles, then the on-shell residue may be seen to be IR-divergent 



IR 



ie 2 f d A k q^p u g^ v 



[p 2 — m 2 ] [q 2 — m 2 ] J (27r) 4 (q ■ k + ie)(p ■ k — ie) k 2 + ie 



(73) 



where we have dropped higher powers of k as they do not lead to IR- divergences (this 
means that we are henceforth restricting ourselves to loop momenta smaller than some 
cutoff). The easiest way to calculate these divergences is to perform the k integral 
using Cauchy's theorem. The four poles are: 

k = ±\k\ =F ie , k = q ■ k/q - ie , k = p ■ k/p + ie . (74) 

The results of the two different sets of poles have different physical interpretations. To 
see this distinction, let us take them one at a time. 

The contribution from the k 2 poles to (^) is easily found to be 

1 1 log A + K el | \ rdW 



(p 2 — m 2 )(q 2 — m?) 8tt 2 \v re \\ \1 — \ v rc\\ J J \k 

where v re \ is the relative velocity between the two charged particles ||. In (f75|) there 
is a logarithmic divergence in the small |fe| limit, i.e., there is a divergent contribution 
from soft virtual photons. This is called a soft divergence. 

The other poles are often neglected as they correspond to divergences in the (unob- 
servable) phase of the Green's function, which cancel in cross-sections. These structures 
yield in flt|) 

-e 2 1 r 1 du f d\k\ , . 

/ / I I (76) 



(p 2 — m 2 )(q 2 — m 2 ) An 2 J-i \v TC \\u — ie J \k 



20 



An IR-divergence is again visible. To perform the w-integral we can now employ the 
relation 

= PV- + m5{u) , (77) 
u — le u 

where PV denotes the principle value. Only the last term in (|77|) contributes as the 
other leads to an odd integral in u. We so obtain 

-ie 2 1 1 rdlkl , s 



{p 2 — m 2 ){q 2 — m 2 ) Aw \v ie \\ J \k\ 

The extra factor of i here, compared to ([75Q , betrays the fact that this divergence 
occurs in the (unobservable) phase. Such singularities are called phase divergences. 
(Note that this is why we chose a pair creation process: if we had taken a scattering 
process then both of the poles in p i ■ k would be in the same half-plane and could be 
neglected.) Now we want to demonstrate that the use of the correct, gauge invariant 
fields removes both of these divergences. 



4.1.2 Perturbation Theory with Dressed Fields 

Using dressed fields means calculating Green's functions of the fields given in fl56|). 
Since the dressings explicitly depend on the coupling, e, we must also expand the 
dressings as well as including the usual interaction vertices. Thus we introduce new 
vertices and hence new diagrams. The diagrammatic rules for dressed Green's functions 
are then just the standard ones augmented by the extra vertices which come from 
expanding the dressings in powers of the coupling. The two factors in the dressing each 
yield a different vertex structure. The Feynman rules corresponding to the dressings 
are given in Fig. 4.3. 



fJL 




Fig. 4.3: The Feynman rules from expanding the dressings. The first vertex comes 
from the soft (x) part of the dressing, and the latter corresponds to the phase (K) 
term. 

The dressings are, of course, dependent upon the momentum of the particle being 
dressed and so here we have defined 

v? = (,+.)"(, -.)•*-*", K--= (n+v) '^ ( ; + % +vYk2 < 79 > 
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where v = (0, v) is the velocity of a particle with momentum p = 7717(1, v). 

Since the dressed fields are, by construction, gauge invariant, so are their (con- 
nected) Green's functions. In the connected vertex function, there are two further 
diagrams with extra interaction vertices, Fig.'s 4.4.b-c. The remaining diagrams of 
Fig. 4.4 come from expanding both parts of the dressing. 




(d) ' (e) (/) 




(9) (h) (i) 




Fig. 4.4: All one-loop Feynman diagrams in the pair creation process which contain 
IR-divergences. Diagrams (a) — (c) are covariant; (d) — (/) involve the perturbative 
expansion of the x term in the dressing; (g) — (i) comes from expanding the K term; 
finally the diagrams (j) and (k) are cross terms from expanding both dressing struc- 
tures. 

Our procedure is to extract the IR-divergences in the on-shell residue after we have 
taken out a pole for each external leg. For simplicity we work in Feynman gauge, 
however, the gauge invariance of our dressed Green's functions will be fully apparent 
in our final results. We will employ dimensional regularisation since it preserves gauge 
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invariance. We may neglect the sorts of diagram shown in Fig. 4.5: those of type (a) are 
massless tadpoles (and so vanish in dimensional regularisation), type (b) do not yield 
IR-divergences when we extract the poles, which reflects the spin independence of the 
infra-red structure, and diagrams of type (c) are also found not to yield IR-divergences 
in the residue. 



Fig. 4.5: Classes of diagrams which do not contribute IR-divergences to the residue. 
The square vertex here and below signifies that the generic contributions of both parts 
of the dressing are meant. 

Let us now consider the contribution from the x P ar t of the dressing. Since the other 
term in the dressing is itself gauge invariant, this part, taken together with the covariant 
diagrams, must be gauge invariant. The diagrams 4.4.e-f have two poles already, but 
Fig. 4.4.d appears not to have any poles: it is, however, already IR-divergent even 
off-shell. This has nothing to do with the usual IR-divergences which, we recall, arise 
when we go on-shell. We can, however, extract poles from such diagrams as we now 
describe. 



4.1.3 Factorisation 

This sort of off-shell IR-divergence appears when one or more photons are exchanged 
between two dressings. We will call such diagrams rainbow diagrams. When we 
extract poles associated with the external legs, we find the residue has an on-shell 
IR-divergence. Together with the other diagrams of Fig. 4.4, the divergences from 
the rainbow diagrams will then (if our predictions of the IR-finiteness of the physical 
Green's functions is correct) cancel the soft and phase divergences of the covariant 
(dressing independent) diagrams. This procedure makes much use of the following 
algebraic identity: 
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(a) 



(b) 




1 1 



1 + 



2p ■ k - k 2 



(80) 



{p — k) 2 — m 2 p 2 — m 2 



{p — k) 2 — m 2 
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The integrand from the Feynman rules for Fig. 4.4.d may be thus rewritten as (we do 
not write superfluous factors) 



/ 



d 4 k 



v P -v q 



(2tt) 4 V p ■ kV q ■ kk 2 [(p + kf -m 2 + ie] [(q - k) 2 

d 4 k Vn-Vn 1 1 



(27r) 4 Vp ■ kV q ■ k k 2 p 2 — m 2 (q — k) 2 — m 2 + ie 



m 2 + ie\ 



1 - 



2p ■ k + k 2 



(p + k) 2 — m 2 + ie 



The first term in the square bracket here is the only relevant one. Before we consider 
it, let us show that the second term does not contribute to the residue: extracting the 
pole in l/(q 2 — m 2 ) we obtain from this term 



d 4 k 



v P -v q 



k 2 - 2p ■ k 



1 + 



2q ■ k - k 2 



(q — k) 2 — m 2 + ie 



[p 2 — m 2 ] [q 2 — m 2 ] J (27r) 4 V p ■ k V q ■ k k 2 (p + k) 2 — m 2 + ie 

(82) 

We have extracted the two poles and may now go on-shell: but the square bracket 
here, and hence the contribution of this second term, vanishes on-shell. Returning now 
to ([81]) the contribution of the first term is easily found to be 



1 



d 4 k 



V p -V q 



[p 2 - m 2 } [q 2 - m 2 } J (2tt) 4 V p ■ kV q ■ kk 



1 + 



2q ■ k - k 2 



(q — k) 2 — m 2 + ie 



(83) 



However, the one in the square bracket here is just a massless tadpole and thus vanishes. 
On-shell the second term becomes —1 and we obtain the final result that, as far as IR- 
divergent terms in the residue are concerned, we have 



/ 



d 4 k 



v P -v q 



1 



(2tt) 4 V p ■ kV q ■ kk 2 [(p + k) 2 - m 2 + ie] [(q - k) 2 - m 2 + ie] 



x 



d 4 k 



V p -V q 



(2Tx) 4 Vp-kV q -kk 2 



\p 2 — m?\ [q 2 — m 2 ] 

i.e., the rainbow line has been stripped off the original integral to yield a factor of 

d 4 k 



J4) 



a 



pq 



v P -v q 



(2ir) 4 V p -kV q -kk 2 



^5) 



times the Feynman rules for the diagram without this line. It may be demonstrated 
that this property of the factorisation of rainbow lines is completely general. For a 
diagram with n rainbow lines, as far as the soft divergences in the on-shell residue of 
the poles are concerned, the rainbow lines may be stripped off and replaced by (C pq ) n . 
A diagrammatic proof of the factorisation property will be published elsewhere. We 
will make extensive use of this property in the all orders proof below. 
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4.1.4 Tadpoles 



With the factorisation property we are able to calculate all of the one-loop diagrams. 
However, it seems worthwhile to first remark on the identification of IR- and UV- 
singularities since, as is well known, in dimensional regularisation we can exchange IR- 
and UV-divergences because massless tadpoles vanish in this scheme. We have used 
this cancellation above in fl83|). Using this identity the badly defined, soft, off-shell 
divergences were removed and replaced by terms which yield IR divergences only when 
we went on shell. These on-shell singularities, as we will soon show, help to cancel 
the soft on-shell divergences of the other diagrams. This replacement of ill-defined, 
off-shell singularities in rainbow diagrams is the only place in our calculations where 
we drop tadpoles. At no stage do we alter structures which first develop soft or phase 
divergences only on-shell. 

However, there are other massless tadpole diagrams associated with dressed Green's 



functions, as well as the tadpoles which we dropped above in Eq. |83|. Let us now 
consider their interplay — firstly for the dressed propagator. Fig. 4.6 shows the generic 
diagrams which enter here: 



(a) (b) (c) (d) 

Fig. 4.6: The classes of diagrams in the one-loop dressed propagator. 

Here there are massless tadpole diagrams, Fig. 4.6.d, and rainbow diagrams, Fig. 4.6.c. 
The propagator has, of course, no phase divergence, so let us solely consider the con- 
tribution of the soft part of the dressing hereQ. Using fl30|) to strip off the rainbow line 
we obtain from this diagram (for the soft part of the dressing) the massless tadpole 
integral 

^ f d 4 k V p 2 -i 
p 2 -m 2 J (2vr) 4 (V p -k) 2 k 2 ' K ' 

which we can drop in dimensional regularisation. However, the contribution of the 
two massless tadpole diagrams, of the type shown in Fig. 4.6.d, is rapidly found from 
the Feynman rules to be equal to ( p6|) but with the opposite sign and so cancels it 

n It is easy to check that the phase part of the dressing does not bring in any soft divergences 
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exactly^. The remaining integrals conspire to yield an infra-red finite propagator as 
has been shown elsewhere in detail GO, f4l . 



This exact cancellation does not carry through to higher Green's functions. In 
general a gauge invariant subset of tadpoles survives. This appears to be a subtlety 
concerned with taking products of distributions. It has been suggested that such 
problems may perhaps be removed by using the Hertz formulation of QED HT^| . The 
removal of these tadpoles in a general fashion is a worthy object for further study, here 
we content ourselves with noting that they vanish in dimensional regularisation and 
that we do not drop tadpoles connected with on-shell infra-red divergences. 



4.1.5 Cancellation of the Soft and Phase Divergences 

Putting together the one- loop diagrams relevant to the soft structure, i.e., Fig.'s 4.4.a-f, 
extracting a pole for each of the external legs and going on-shell, we obtain the following 
terms in the residue which are by power counting IR-divergent: 



<m = 

J p I \ p ■ K V q • K I K" \Vp-K (/ • K I 

(87) 

where we have removed a factor of ie 2 /[(p 2 — m 2 )(q 2 — m 2 )]. This is a gauge invariant 
result: replacing the Feynman gauge propagator by any more general form necessarily 
brings in a factor of or k v which vanishes on multiplication into the brackets in (|S7j). 
Similar gauge invariant structures are obtained in the study of the dressed propagator 
and in other physical vertex functions. We now have to show that the soft divergences 
seemingly apparent in (pTT) in fact cancel. 




The cancellation may be demonstrated by performing the integrals explicitly, using 
the methods outlined above. However, it is simpler to realise that 

V£ _ {y + vy(ri-v)-k-k» 



Vp-k k 2 - (fe • v) 2 

which, at the soft pole, k 2 = 0, we may write as 

Vp (t] + vY(t]-v) -k- HP 
V p -k ~ (k ■ r]) 2 - (k ■ v) 2 



i9) 



12 Note that the massless tadpole diagrams of Fig. 4.5.b do not cancel, but they, of course, do not 
alter the on-shell dependence of the IR-structure 
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We may drop the k^ term in the numerator here (cf, the argument for the gauge 
invariance of Eq. [37|). So effectively we have 



V£ (rj + vYij] - v) ■ k (tj + vY p 



V^> • (k ■ rj) 2 — (k ■ v) 2 (rj + v) ■ k p ■ k ' 

where it is important to note that we have taken p to be at the correct point on the 
mass shell, p = 7727(77 + v ), which defines the appropriate dressing. (If we were to go 
on- shell at a different point, this equality would not hold and the soft divergences do 
not cancel.) 

In all cases of gauge invariant Green's functions, we are thus able to replace 

" ' -d T7 ^- = ^ T , (91) 



V p ■ k p ■ k V g ■ k q ■ k 

as far as the soft divergences are concerned. This equality makes the cancellation of 
the soft divergences in (8l) and other dressed Green's functions immediately apparent. 



Up till now we have not included the other gauge invariant part of the dressing. It 
is straightforward to check that the contribution of this structure to the pair creation 
process is gauge invariant. It may also be demonstrated that it does not spoil the above 
cancellation of soft divergences. (Some of the individual diagrams which involve this 
part of the dressing, Fig.'s 4.4.g-k, contain soft divergences, but overall they cancel.) 
The role of this structure is, we expect, to cancel the phase divergence of Fig. 4. 4. a. The 
only phase divergent contribution from the entire dressing comes from a single diagram 
involving only the -fT-structure: Fig. 4.4.g. From the Feynman rules we obtain from 
this figure after factorisation 



d 4 k Wl 9llv f d A k ( V + V ).( v + V ')k 2 



J (2tt) 4 V p - kV q -k k 2 J (2vr) 4 V p ■ kV q ■ k[(r, + v) ■ k - ie][(r, + v') ■ k + ie] ' 

where we have not written the two poles — ie 2 /([p 2 — m 2 ][q 2 — m 2 ]) and have thrown 
away terms which do not yield phase divergences. This result may be readily seen to 
cancel exactly with (^) when we go on-shell at the right point. 

To summarise the one-loop results: the dressed Green's functions are, by construc- 
tion, gauge invariant. If we go on-shell at the points on the mass shells such that the 
momenta of the external lines and the velocities that define the various dressings agree, 
then the soft and phase divergences cancel. We now want to show that this holds at 
all orders. 
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4.2 All Orders 



The physical dressed propagator was shown elsewhere f|l| to be IR-finite at one 



loop. However, we can show that this holds at all orders through a slight extension of 



the work of Jackiw and Soloviev | 43| . 



These authors used low energy theorems and spectral representations to study the 
behaviour of the propagator around the pole, p 2 = m 2 , in various gauges. In covariant 
gauges they regained the usual result (|15D , while in non-covariant gauges there is a 
pole if the following integral vanishes (this is their Eq. 3.41) 



" 2 / ke* kk 9(k )6(k 2 )^Il% , (93) 



1 R 



(2nfJ v uy v '(Hfe) 



here r corresponds to the mometum of the matter field and II is the photon propagator 
up to a factor of i/k 2 . Now our dressed propagator, for a charge with momentum 
r = m7(r/, v) corresponds to the usual matter propagator in the gauge where 

n K-- g + k^Tv) (H v + v)) 2 \- K +2Vr ' k \ ■ (94) 

taking the S(k 2 ) factor of the integrand into account, it is easy to see that at the right 
point on the mass shell ( |9"3] ) vanishes and the propagator has a pole. This was pointed 
out explicitly for the static case (Coulomb gauge) in Ref. [[43[ . 



4.2.1 The Dressed Propagator 

Since we now know that the dressed propagator has a pole, we see that the wave 
function renormalisation constants are IR-finite. We can therefore multiply the dressed 
Green's functions by them without introducing new IR-divergences. We will use this 
freedom to show the cancellation of soft and phase singularities in dressed Green's 
functions at all orders. Our approach is diagrammatic. First we shall consider the 
types of diagrams which can introduce IR-divergences in the propagator, in this way 
we can describe the structure of the wave function renormalisation constants. 

The non-dressed propagator, i.e., just the covariant diagrams, can be represented 
at all orders in the following way: 



p2 _ m 2 
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When we take dressings into account, these chains must be supplemented by various 
possible structures which we now enumerate. The first is 



2i£ 



P 



where the 2 accounts for the possibility that dressing corrections, S, can be indistin- 
guishably attached at either end in the scalar theory. 

There may of course also be such dressing corrections at both ends. This is dia- 
grammatically 



P 



77T 



Taking all of these possibilities together we see that these contribution to the prop 
agator have the form 



p 2 — m 2 — S 



(95) 



Of course we can also exchange one or more photons from one dressing to the next. 
This can be done for each and every one of the above sets of diagrams. We so obtain 
the following rainbow diagrams 





Using the factorisation property, summing over all possible numbers of rainbow lines, 
and taking into account the 1/n! symmetry factor which accompanies a diagram with 
n rainbow lines, we can re-express this algebraically as 



p- 



where — C pg was defined in Eq. EB 



77T 



• exp {-C pp ) , 



(96) 
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Of course there can be other diagrams such as 




where the dressing corrections overlap in a non-rainbow form, i.e., diagrams which are 
not one-particle irreducible even after rainbow lines are stripped off. Such diagrams 
do not yield poles and so we can drop them here. We neglect loops of matter fields 
since these also remove the soft divergences. Eq. ^ is therefore our final result for the 
infra-red divergences in and the pole structure of the physical propagator. 

We see that the dressings factorise into a covariant and a non-covariant part. For the 
full non-covariant wave function renormalisation constants, we thus have the attractive 
result that 

Zf = Z 2 cov (l + S) 2 exp(-C pp ). (97) 

Although only a partial exponentiation of the dressing contribution is apparent, we 
know that it is in fact IR-finite at all orders. However, fl97l) will suffice to show that 
the pair creation vertex is finite to all orders. 



4.2.2 All Orders Pair Creation 

The general class of diagrams with double poles and possible IR-divergences in the pair 
creation process is 




\J Z2Z2 



where we have chosen to multiply the diagrams by an (IR-finite) factor: the inverse 
square root of the appropriate wave function renormalisation constant (as given by 
Eq. |9?|) for each of the legs. This is useful as it makes the exponentiation of the 
infra-red divergences apparent. 

The diagrams we have retained are those which can yield IR-divergences, i.e., where 
there are possible covariant vertex corrections from one leg to another, covariant cor- 
rections on the external legs, dressing corrections at the end of the external legs and 
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possible rainbow corrections from one dressing to the other. The use of a black blob, 
• , here denotes that there may or may not be a dressing correction (E) at the ends of 
one or both of the lines. 



Other diagrams (e.g, where a line from a dressing connects to a covariant interaction 
vertex on the other line) will not give two poles or will be infra-red finite. 

We now factorise the rainbow dressings, this means (summing over all possible 
rainbow lines and including the 1/n! symmetry factor for a diagram with n such lines) 
that we may write 




Since there may or may not be a dressing correction at the end of these lines, we 
may write these end factors as (1 + E). They clearly then cancel the E dependence 
from each of the external leg wave function renormalisation factors as given by fl97[). 
Diagrammatically we thus have 




and we see that the dressing effects have exponentiated. Since the covariant diagrams 
which are left over are well known to exponentiate, we see that all the soft and phase 
effects exponentiate in the residue of the double poles. We have seen that they cancel 
at one loop, this now trivially implies that the dressed process is IR-finite at all orders 
at the level of Green's functions. We stress that this holds for both soft and phase 
divergences. The extension of this approach to scattering and to higher vertices is 
straightforward. 

We have thus seen that the dressings we have introduced in Sect. 2 remove soft and 
phase divergences in massive Electrodynamics already at the level of Green's functions 
and permit a particle interpretation. The immediate question now is whether these 
perturbative successes can be extended to QCD. In particular, we need to clarify the 
generalisation to a different class of mass singularities, collinear divergences, which 
arise in QCD. Since collinear divergences also characterise QED with massless charged 
particles, the next section is dedicated to a study of that theory. 
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5 Massless Charges 



Massless charges are an essential ingredient of the standard model. In QCD the colour 
charge of the massless gluon results in its self-interaction which in turn leads to asymp- 
totic freedom. In addition, gluonic bremsstrahlung plays a dominant role in jet forma- 



tion [P]. However, such massless charges are poorly understood, even at the level of 
cross-sections, since the Block-Nordsieck procedure of summing over degenerate final 
states fails to yield finite results [f45|, ^6[. Massless charges can be modelled in QED 
by taking the limit of vanishing electron mass. As one would expect, QED then has 
many new and unusual features. In particular, massless electron pairs can now be 
created by expending arbitrarily small amounts of energy. This freedom leads to a 
new type of singularity, the collinear divergences. In perturbation theory these express 
themselves through the divergence of on-shell massless propagators. The inverse prop- 
agator, (p — k) 2 — m 2 , vanishes in the m — > limit if p is on-shell and k is parallel to 
p. (Note that k is not necessarily small.) For further details see ||. The matter pairs, 
which can now be easily created, screen the initial charge and the effective coupling 
vanishes [0. Here we will see that the argument of Sect. 2 for the vanishing of the 
asymptotic interaction Hamiltonian, if correctly dressed matter is used, also applies in 
the massless limit. Finally we discuss some aspects of the construction of dressings in 
this limit. 



5.1 Collinear Asymptotic Dynamics 

The asymptotic behaviour of the interaction Hamiltonian for massive QED was dis- 
cussed in Sect. 2. In the massless case we will show that the interaction Hamiltonian 
in the distant past or future has a far richer structure than in the massive case. 

The interaction Hamiltonian, 7ii nt (t), is given by (|8|), where the conserved matter 
current and the free field expansions for matter and gauge fields are still those of (§) 
and (|i0|). The all-important difference is that, since m — 0, the energy in the massless 
case is given by E p = \p\. 

When the expansions (§) and ( fLOf) are substituted into (H), then, of the eight possible 
terms, the two involving a^vj^u and ajjl^v will have a time dependence of the form 
e tat with a being either positive or negative. These cannot survive for large values 
of t and so do not appear in the asymptotic limit. This leaves six terms which are 
potential survivors, of which two vanish in the massive case. A closer look at some of 
these structures will show how to evaluate all the terms. 

The first term to be considered in detail is one of the "off-diagonal" (in the spinors) 
terms: 

d?xd 3 kd?pd 3 q a„ (fc) fet ^ r ^^ v r^ e -ik. Xe i q .* e ip.* (9g) 
(2tt) j 2u k J4E p E q 
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Integrating out the cc- integral will give a term involving 5(q + p — k). Following this 
by integrating out the p integral, the expression ( |98|) becomes 

d*kd3q a jk) ftf f _ ^ _ y^^,^ (gg) 

(2tt) 6 2u k J4E q E k - q 



If this is to survive as £ — > oo then the coefficient in the exponent must vanish. This 
is equivalent to demanding that uj k = E q + E k - qi or \k\ — \q\ + \k — q\. Since these 
three vectors represent the three sides of a triangle, plainly a solution to this equation 
will exist if and only if the vectors k and q are parallel and \k\ > \q\. In this region, 
where the photon is collinear with the matter field, this contribution to the asymptotic 
interaction Hamiltonian does not vanish. The momentum, k, is no longer restricted to 
the soft region, k = 0, but must be larger than the matter momenta, since this term 
corresponds to pair creation. 

The second term to be examined is one of the "diagonal" terms and has the form 

- e / d3xd ^ d 'P d ^ a f) d(g> 8)d \fa r)^)/^^ . (100) 
J (2vr) 9 2uj k ^JAE p E q 

Integrating out the x again gives a delta function, this time S(p — k — q). Integrating 
out the q integral now yields 



d 3 k d 3 p 



(2tt) 6 2ujAE p E p 



^ffl d{p - k, s)S{p, r)v s {p - k^v'ip)^-^-"-^ . (101) 



The exponent in ( |101| ) must vanish for large t and this implies that E p — E p _ k — uj k 
must vanish. This is similar to the previous case except that now we must have 
\p\ = \k\ + \p — k\. Reasoning as before, we find that p must be parallel to k but now 
\p\ > \k\, i.e., E p > u>k- This clearly corresponds to photon production. 

If the six terms are evaluated using (|128[) then the final form of the asymptotic 



interaction Hamiltonian in the massless case is |48 

/d k 
— [ afl (k)J^(k,t) +h.c], (102) 

where 'h.c.' denotes the hermitian conjugate of the first term in brackets, and 

J&(k,t) = J ^^(Pscattfefc) +p prod (p, k)) ■ (103) 

The two structures in the asymptotic current are respectively 

Pscattfo k) = YH^iP, S)b{p -k, S )~ d\p, S)d{p - k, s )yt[Ep-E p _ n) ^ 
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in the region Uk < E p , and 

P P rod(p, fe) = E fot (P» s ) rft ( fc - P» r )^ t<T • g£V*^ +E ^-<^ (105) 



in the region u>k > E p . The two terms in ( |102|) have different physical interpretations: 
the p SC att part corresponds to photon radiation and includes the soft region which was 
responsible for the non-vanishing of the usual soft asymptotic interaction Hamilto- 
nian, (|TTD . However, the photon momenta is only required to be collinear and not 
extremely soft. The other term, p pro d, is completely new and corresponds to the pro- 
duction of massless matter pairs with momenta less than that of the initial photon. 
Even though photons cannot radiate other photons, these two structures in this model 
theory already show the basic processes underlying the collinear production of gluons 
and quark-antiquark pairs in jet creation. 



We now note that from the form of ( 102 ), and in particular the explicit factor 



in the asymptotic current, the solution to the dressing equation, (|24]) or ( p5f) , will, also 
here in the massless theory, correspond to a particle since the asymptotic interaction 
Hamiltonian vanishes at the correct point on the mass shell. We thus predict that the 
Green's functions of the dressed fields will also be free of collinear divergences. 

In the context of this discussion of the asymptotic interaction Hamiltonian, it should 
be noted that, for massive electrons, if we allow the photon a small mass there is no 
momentum configuration such that the exponential survives at asymptotic times. In 
perturbation theory it may be easily seen that a small photon mass regulates the infra- 
red singularities. However, if the mass of the photon becomes large enough to open 
decay channels into matter pairs, then once again the interaction picture breaks down 
and this results in apparently gauge dependent S-matrix elements (see, e.g., [§j|). 

We now turn to the construction of the solutions of the dressing equation in the 
massless limit. 



5.2 Collinear Dressings 

In the massless theory, the dressing must still fulfill the requirements of gauge invari- 
ance (^) and the dressing equation (25), where now (77 + v) 2 = 0. Following the 
procedure outlined in Sect. 2 for solving these equations, we will require the commuta- 



tors of the potential in the theory described by the asymptotic current, (|1(J3[) . As we 



saw for the massive case, the potential can now be written as in fl4"7). However, in this 
massless theory, the asymptotic currents no longer commute and the commutators of 
the theory are no longer those of the free theory. This makes the construction of the 
dressing of massless charges somewhat more difficult. Rather than discuss the details 
of how to extend the analysis of Sect. 2, we will now examine the v — > 1 limit of the 
massive case and then study, in this limit, the route taken by Dirac to the dressings. 
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The naive massless limit of a dressed charge ( p)B"D moving in the x 1 -direction results 
in integrals of the form 

1 f „ (d 2 A 2 + d 3 A 3 -Ei)(x°,z) , , 

An J \x x - z x \ 

Clearly this is ill-defined, and although l/lx 1 — z 1 \ can be defined in terms of generalised 
functions , the serious problems with this naive limit cannot be circumvented in this 
fashion. Further evidence that a naive approach to the massless theory will not suffice, 



comes from the detailed perturbative calculation of the one-loop propagator (40[ . 



The wave-function renormalisation constants found in these papers diverge as v — > 1. 

To obtain a deeper insight into the form of the dressings, we can follow Dirac's lead 
and construct a dressing that reproduces the electric and magnetic fields for a massless 
charge. The electric field of an ultra-relativistic charge contracts in the direction of 
motion. For a massless charge this contraction leads to a singular field configuration. 
It has been argued that, for motion in the x 1 -direction, the electromagnetic fields 
are 0, ^ 

— ex l 5(x° — x 1 ) 
Ei = 0, E l = - 



B 1 = 0, B. 



2vrr 2 

(107) 

e lj x j 5(x° — x 1 ) 



2nr 



2 



where i, j = 2, 3 and r 2 = (x 2 ) 2 + (x 3 ) 2 = \x±\ 2 . This result should follow from the soft 
dressing term. 

We argue that this part of the dressing should have the form 

X i(x°, x) = ^J d 2 z ± (d 2 A 2 + d 3 A 3 - E^x , x\ z x ) log \x L - z L \ + g(x) (108) 

where g(x) is a function which is harmonic in x ± , i.e., g(x) is in the kernel of the two 
dimensional Laplacian. Although this is not the naive limit (|106|) , it can be understood 
as arising from that part of the dressing equation which describes the soft dynamics, 
Q ' 9(x) — G • A. For this configuration, the equation becomes 

((1 - v 2 )d\ + d\ + dl) Xv {x\ x) = ((1 - v^A, + d 2 A 2 + d 3 A 3 - vE^x , x) . (109) 

Now in the v — > 1 limit, one sees that the limiting value of Xv, which we write Xi, must 
satisfy 

(dl + d 2 3 )xi(x°, x) = {d 2 A 2 + d 3 A 3 - E{){x\ x) , (110) 

with solution ( |108| ). It may be checked that the dressing constructed out of this soft 
term, with g(x) set to zero, indeed yields the electric and magnetic fields, ( |107| ). This 
result strengthens our claim that a careful analysis of the dressing equation will allow 
us to construct massless charges. 
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6 Discussion 



In this review we have demonstrated that it is indeed possible to construct relativistic 
charged particles. Let us now recall what we have seen. 

Our starting point was the non-vanishing of the asymptotic interaction Hamilto- 
nian which characterises gauge theories. The misidentification of the asymptotic inter- 
actions of the free Hamiltonian with the true asymptotic dynamics of gauge theories 
causes the infra-red problem. Since the interaction Hamiltonian does not vanish at 
large times, Gauss' law tells us that charged particles are not just the matter fields 
of the Lagrangian: physical particles like the electron are always accompanied by an 
electromagnetic dressing. In practical calculations this expresses itself in the lack of a 
pole structure in the on-shell Green's functions of the Lagrangian matter fields and in 
divergences in S-matrix elements. 

Any physical degree of freedom must be gauge invariant. For charged particles this 
means that the dressing together with the matter field at its core must be locally gauge 



invariant — this translates into an equation for the dressing, fllq) . This requirement, 
though, does not suffice to construct charged particles: any gauge invariant solution 
is in principle a physical degree of freedom, but it is not necessarily one that physics 
chooses to use. To restrict ourselves to the solutions which are physically relevant, we 



require a second equation, fl25[) . This latter relation was deduced from demanding that 
the velocity of an incoming or outgoing charged particle is well defined and the dressing 
must take this into account. The form of the asymptotic interaction Hamiltonian is such 
that it vanishes in the propagator of a correctly dressed particle, i.e., one which satisfies 
(|T6| ) and (25). It follows that dressed charges are then described by the free Hamiltonian 



and a relativistic particle description is indeed possible. In Sect. 2 dressings which solve 
these equations were explicitly constructed. These are physical degrees of freedom with 
a specific physical interpretation. The dressings factorise into two parts, each of which 
was interpreted as playing a different role in the infra-red physics of QED. Sect. 3 then 
illuminated the role of velocity in any description of charged particles and showed how 
the form of the dressing could be derived from the theory of heavy charges. 

These arguments have been tested by explicit perturbative calculations described 
in Sect. 4. The interpretations of the two different terms in the dressing were such 
that they were expected to each cancel a different type of infra-red divergence. The 
constructions passed these one-loop tests with flying colours and a gauge invariant 
pole structure was obtained. This was then generalised to an all orders proof of the 
cancellation of infra-red divergences in the on-shell Green's functions of QED with 
dressed fields. These results give us great confidence in the requirements (|16|) and (p5|), 
in the solutions Q5"BD and also in the interpretation we associate with these physical 
degrees of freedom. 

Finally in Sect. 5 collinear divergences were studied in the framework of QED with 
massless fermions. Here the structure of the interaction Hamiltonian which survives at 



36 



large times is much richer than in QED with massive matter. However, it still vanishes 
if correctly dressed matter, satisfying fl25|), is used. This implies that the on-shell 
Green's functions of the solutions to (|T6|) and fl25|) in the massless theory will also be 
free of all infra-red singularities, including collinear divergences. The construction of 
these dressings was considered. 



6.1 Charges in the Standard Model 

Most of this paper has been given over to QED. The other interactions in the standard 
model of particle physics are also described by gauge theories, however, the charged 
particles in these theories are very different to those of the unbroken abelian theory. In 
this subsection we will sketch some of the most important differences which characterise 
the charged particles of the weak and strong nuclear forces. 

The asymptotic Hamiltonian in QCD also does not reduce to the free one. However, 
there is a new problem with defining charges in non-abelian gauge theories such as 
QCD. The colour charge 

Q a = [ d 3 x(JS(x) - tf c E<>(x)A*(x)) , (111) 



in sharp contrast to that of QED, is not invariant under gauge transformations. It 
is therefore natural to wonder how we may speak of colour charged particles. But on 
physical states, where Gauss' law holds, the charge may be written as 

Q a = - [ d 3 xdiE?(x). (112) 

Using Gauss' theorem it follows that the colour charge expressed in this way is in fact 
invariant under gauge transformations which reduce, in a directionally independent 
manner, to elements of the centre of the group at spatial infinity. Thus the concept 
of colour charge is only meaningful for locally gauge invariant fields and then only if 
this restriction on the allowed gauge transformations is imposed. Since the Lagrangian 
matter fields are not gauge invariant, coloured quarks in QCD are necessarily dressed 
by glue HI. 



The dressing equation ( ff5| ) of QED may be directly generalised [[^J to QCD. The 
non-abelian solutions will be asymptotically described by a free Hamiltonian provided 
that the field transformations required are, in fact, admissible. To the extent that they 
may be constructed, it follows that the various mass singularities would cancel at the 
correct points in the mass shell of the dressed quark fields. 

However, there is a fundamental obstruction to the construction of dressed quarks 
(and gluons) which follows from the nature of non-abelian gauge transformations. Any 
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gauge invariant description of a quark could be used to construct a gauge fixing. How- 
ever, the boundary conditions which must be imposed on the allowed gauge transfor- 
mations in order for colour to be associated with the physical degrees of freedom, are 
such that the Gribov ambiguity |53|, EM holds. There is then a fundamental, non- 



perturbative limit on the construction of gauge-invariant, coloured charges. Thus the 
true degrees of freedom in QCD outside of the perturbative domain do not include quarks 
and gluons. Colourless gauge-invariant fields can of course be constructed, but quarks 
and gluons are confined. 

The picture of confinement which emerges from studying the construction of coloured 
charges is as follows. When a QQ-system is separating, but the matter fields are still 
at a short distance from each other, the interaction potential is essentially Coulombic 
and the dynamics is described by perturbative dynamics. The short-distance, Coulom- 
bic inter-quark potential may be described using low-order perturbative solutions to 
the dressing equations for each of the individual quarks. Deviations from this such 
that a confining potential arises in a QQ-system are expected to come from a 'mesonic 
dressing' which does not factorise into two parts. 

Finally, it might be objected that the non-abelian theories which underly the weak 
interaction do not lead to confinement of weakly charged particles such as the W or 
indeed the electron. This is easily understood: in spontaneously broken gauge theories 
we are entitled to use the Higgs sector to dress charged particles |55[. (The ability to 
choose the unitary gauge circumvents the gauge fixing ambiguity.) In this way gauge 
invariant solutions corresponding to weakly charged particles may be constructed. 



6.2 Open Questions 

This review has sketched out a systematic approach to the construction of charged 
particles. The qualitative results to date of the program are fully in accord with phe- 
nomenology and the calculational tests to which the charged particles have been sub- 
jected have all supported the methods. But there are still a large number of unanswered 
questions ranging from calculational procedures to the extension of the applicability of 
this programme to physics at finite temperature and density. We will now conclude by 
listing some of the most pressing tasks. 

a) The above proof of the infra-red finiteness of the dressed on-shell Green's functions 
needs to be extended to full calculations of these Green's functions, in particular the 
UV renormalisation of the n-point functions of these composite operators must be 
performed. The dressed propagators of both fermionic and scalar QED have been 
carried out at one-loop. This now needs to be extended to higher loops and to vertex 
functions. 

b) The perturbative tests need to be extended to collinear divergences. Massless QED 
is the natural testing ground here. First the dressing equation needs to be solved within 
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the framework of the remarks of Sect. 5. Then the solutions of the equivalent dressing 
equation for QCD should be constructed and tested. 

c) Although a brute force, direct perturbative solution of the dressing equations for 
the non-abelian theory is feasible at low orders, a systematic and practical approach 
to the construction of a gauge invariant dressed quark fulfilling the dressing equation 
of QCD for an arbitrary velocity is needed. In the appendix to Ref. || an algorithm 
was presented with which the calculation of a particular gauge invariant dressed quark 
solution to any order in perturbation theory became rather simple and some of the 
first terms in this perturbative construction were given. In this context we also refer to 



Ref. []5(| where all orders expressions apparently corresponding to the gauge invariant 
extension of this term were presented. This corresponded to the gauge invariant ex- 
tension of the soft term in the QED dressing of a static charge. Such work needs to be 
generalised to both terms and indeed to arbitrary velocities. Furthermore the success of 
phenomenological constituent quark models in describing hadronic structure strongly 
indicates that, within the overall constraint of confinement, some non-perturbative in- 
put into the construction of dressed quarks should be possible. Here we are thinking 
especially of the role of instantons and condensates in chiral symmetry breaking. We 
stress that it will not be possible to incorporate all the non-perturbative, topological 
aspects of QCD into dressings of individual quarks or gluons. This, we have argued 
above, is how confinement makes itself manifest. 

d) The importance of the perturbative chromo-electric and chromo-magnetic dressings 
is that they determine the short distance interactions between quarks and have impli- 
cations for jet physics. We would urge, e.g., a comparative study of the distribution of 
glue in the dressings around quarks and gluons which could shed light on the different 
development of such jets. The fruitful concept of parton-hadron duality |57[] could we 
suggest be potentially replaced by a refined version of a duality between (perturba- 
tively) gauge invariant, dressed colour charges and the resulting physical hadrons. 

e) The qualitative proof of quark confinement which was described above needs to 
be quantified. How can the hadronic scale be found in this way? That quarks are 
not part of the physical degrees of freedom of QCD will only become apparent at 
larger distances where that part of the non-perturbative dynamics which is sensitive 
to the Gribov ambiguity becomes first significant and finally dominant. This is, of 
course, a tough non-perturbative calculation. One natural technique is the lattice, 
where there has been some work on studying the gauge fixing problem. It will be 
important here to distinguish between lattice artifacts and the true, physical limitation 



on gauge fixing pq, 59 . We also suggest here that the construction of dressed charges 



in monopole backgrounds be studied. 

f ) In view of the difficulties inherent in non-perturbative calculations, phenomenolog- 
ical modelling of dressings is desirable. Thus we feel that the construction of dressed 
quarks with phenomenologically desirable properties, e.g., running masses should be 
investigated. A corollary of this last point is the construction of mesonic dressings 
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which lead to phenomenological inter-quark potentials. The stability of such model 
dressings should then be tested using variational methods. We recall from the intro- 
duction the instability of the (confined) 'e + e~ '-system where the two matter fields are 



linked by a string [21 



g) One of the most important questions in particle physics is how mass is generated. 
The construction of charges in theories with spontaneous symmetry breaking deserves 
further, quantitative investigation. Both perturbative and non-perturbative effects 
should be studied. A major question closely related to this is how one should describe 
unstable charged particles, such as those occurring in the weak interaction. It is well 
known that there are problems with obtaining gauge invariant results in the presence 
of such fields [|9|, ^l|]. We have seen, at the end of Sect. 5.1, the non- vanishing of 
the asymptotic interaction Hamiltonian in such theories. We urge the development of, 
in some sense dressed, admixtures of fields which would fulfill the following properties: 

i) if the coupling is set by hand to zero it should reduce to the Lagrangian field, 

ii) they must be gauge invariant and iii) the asymptotic interaction Hamiltonian should 
vanish for the propagator of these constructs. 

h) To what extent can we talk about charged particles beyond the standard model? 
The construction of descriptions of charged particles in theories such as technicolour 
and with (broken) supersymmetry could cast light on their experimental signatures. 
Here it is important to see if the gauge invariant, dressed charged particles are still 
related by supersymmetries, or if this just holds for the Lagrangian fields. Furthermore 
we strongly suggest a study of the construction and dynamics of charged particles 
in unbroken supersymmetric models so as to clarify the role of (supersymmetrically) 
charged particles in the Seiberg- Witten description of confinement in such theories |J2] . 

i) Constrained dynamics is the mathematical framework for extracting physical degrees 
of freedom in theories such as QED and QCD. In modern formulations BRST symmetry 
is used to single out gauge invariant, local fields. However, as will have become appar- 
ent above, gauge invariance is not enough: one must still make a clear identification 
between the true degrees of freedom and the observed particles. Quark confinement 
in QCD is merely the most obvious example; even in Quantum Electrodynamics the 
infrared problem, and the associated superselection structures labelling charges with 
different velocities, show that this identification is not direct. The interplay between the 
asymptotic interaction Hamiltonian and the physical observables needs further study. 
It should be noted that the BRST method is a local construction and, as we have 
seen, charged fields are necessarily non-local. In this context we point out a further 
symmetry of QED which was noted in Ref. |53| and which singled out a subset of gauge 
invariant fields. We believe that this is only one representative of a class of symmetries 
with whose aid the true degrees of freedom may be interpreted. Finally, the incorpo- 
ration of nonperturbative effects is the outstanding question in mathematical physics 
today. In constrained dynamics this entails extending the usual reduction procedure to 
theories with superselection sectors characterised by the existence of nontrivial surface 
terms and other global structures. A proper understanding of these effects is essential 
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if constrained dynamics is to be able to discuss nonperturbative physics whether it be 
the quark structure of low energy QCD or string phenomenology 
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A Appendix 

To obtain the asymptotic Hamiltonians of Sect. 2 and Sect. 5 we must evaluate the 
terms u 1 '(p)^^ '(q) , If '(p)^v s '(q) , %T {p)^u s (q) and u r (p)j^v s (q) . The results we will 
obtain, and indeed some generalisations, are quoted in Appendix J of Ref. |H . 



Given an on-shell 4- vector, p, so that p = m , then E p is given by E p = y \p\ 2 + m 2 . 
We introduce the notation 

N P = — and p = N p p. (113) 

hip + m 



Let = (1, 0) and = (0, 1). The Dirac spinors for r = 1, 2 are taken 



as 



"' w =-^u r K0 w-M"*) (u4) 



The u r (p) and v r {p) are, respectively, positive and negative energy solutions to the 
Dirac equation. Let us define 

* = (!£)■ < 115 > 

Then we have 

u r (p) = Av r (p), (116) 

and 

u r {p) =v r (p)-f°A-?° . (117) 
Our gamma matrix convention is 
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where a % are the Pauli matrices. It is now straightforward to show that 

tf(p)YV(g) = v r (p)^v s (q), (119) 
u r {p)Yv s {q) = v'\p)Yu s {q). (120) 

This observation is based upon the easily verifiable identity 

7 °A 7 VA = Y . (121) 

Two other useful observations for the evaluation of the above expressions are the fol- 
lowing. The first is 

u r (p)^u s (q) = ±£f[p)(ff + Yd)u s (q) (122) 

which uses the fact that u r {p) is a positive frequency solution to the Dirac equation; a 
similar identity can be obtained for If (p)^f fl u s (q) . The second is the observation that 

ff + = (p + qf + «(p - q)^ v , (123) 

where o^ v = |[7 M , 7"]. Routine calculations now lead to the following results. 



r (ph°u s ( q ) = ^=e rt <T ■ (p + g)e , 

w r (p)7V(g) = -^=|=£ r t([( CT . p )g + (o- . q) p - tp x q f + (1 - p ■ qr)<7 4 )£ a , 

u r (p) 7 V(g) = -^=f t(l + p ■ q + icr ■ p x q)f , 

vTip^u'iq) = -j±-—^(p + q + icr x (p - q))^ s . 



(124) 



The latter pair of results is sufficient to calculate the asymptotic interaction Hamilto- 
nian ( |TT|) of the massive theory. 

In the m = limit, these results simplify in an attractive way. Recall from the 
discussion of Sect. 5 that the asymptotic interaction Hamiltonian only receives a con- 
tribution from the region where the momenta are collinear. As an illustration, we 
shall examine the first of the terms in (|124j) , in the collinear region where we have 
1^1 = l<7l + l& — q\ ( c -f- the paragraph after Eq. ^). In the massless case we have 
N p = 1/Ep, so that p = p/Ep. It is also clear that p 2 = 1 and p • p = E p . We thus 
obtain for massless charges for this momentum configuration 

v r (qh°u s (k-q) = 1 fV-(q+(fc-q))f. 

'N k - q N q 

E^ q ^ a .(±+*pL\ e . (125) 
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Now in the asymptotic interaction Hamiltonian this spinor combination only occurs 
if the momenta are such that k and q are parallel light-like vectors, and using this in 
(|125|) , we see that 

Tf(g) 7 V(fc-g) = ylE q E k _ q (*q-tT? + ^) ( 126 ) 
= 2^E q E k „ q e ] q-<rC- (127) 

The other terms in (|124j) can be found in a similar manner. We confine ourselves to 
listing them. They yield the simple results: 

lf{q)^u s {k -q)= 2q^E q E k _ q ^q ■ tr£° 
UJ\q)^u s {k-q) = 2p^E q E k _ q 5 rs 

These results may be used to extract the asymptotic interaction Hamiltonian for mass- 
less QED. 
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